Dissipative Landau- Zener transitions of a qubit: bath-specific and universal behavior 



Keiji Saito, 1,2 Martijn Wubs, 3 Sigmund Kohler, 3 Yosuke Kayanuma, 4 and Peter Hanggi 3 

'Department of Physics, Graduate School of Science, University of Tokyo, Tokyo 113-0033, Japan 
Department of Physics, Graduate School of Science, 2 CREST, 
JST, 4-1-8 Honcho Kawaguchi, Saitama, 332-0012, Japan 
3 Institut fur Physik, Universitdt Augsburg, Universitdtsstrafie 1, D-86135 Augsburg, Germany 
^Department of Mathematical Science, Graduate School of Engineering, 
Osaka Prefecture University, Sakai 599-8531, Japan 
(Dated: February 6, 2008) 

We study Landau-Zener transitions in a qubit coupled to a bath at zero temperature. A general 
formula is derived that is applicable to models with a non-degenerate ground state. We calculate 
exact transition probabilities for a qubit coupled to either a bosonic or a spin bath. The nature 
of the baths and the qubit-bath coupling is reflected in the transition probabilities. For diagonal 
coupling, when the bath causes energy fluctuations of the diabatic qubit states but no transitions 
between them, the transition probability coincides with the standard LZ probability of an isolated 
qubit. This result is universal as it does not depend on the specific type of bath. For pure off- 
diagonal coupling, by contrast, the tunneling probability is sensitive to the coupling strength. We 
discuss the relevance of our results for experiments on molecular nanomagnets, in circuit QED, and 
for the fast-pulse readout of superconducting phase qubits. 
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I. INTRODUCTION 

Nonadiabatic transitions at avoided level crossings 
play an essential role in numerous dynamical phenom- 
ena in physics and chemistry. They have been stud- 
ied both theoretically and experimentally in various con- 
texts like spin-flip processes in nano-scale magnets ji'2' 3 -'^ 
molecular collisions, 5 optical systems^^ quantum-dot 
arrays, % Bose-Einstein condensates^ the control of chem- 
ical reactions^ and recently in particular in quantum 
information processing, u i 12 : 13 . u : 15 i 16 : 17 : 18 . 19 . 20 

The "standard" Landau-Zener (LZ) problem describes 
the ideal situation in which the dynamics is restricted to 
two levels that are coupled by a constant tunnel matrix 
element and cross at a constant velocity. The quantity 
of primary interest is the probability that finally the sys- 
tem ends up in the one or the other of the two states. 
This classic problem was solved independently by sev- 
eral authors in 1932j 21 i 22 ' 23 i 24 In quantum devices, not 
only the transition probability but also the nonadiabatic 
relative phase (Stokes phase) between the two states is 
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This phase leads to observable interfer- 



ence effects, for example in superconducting qubits, 17 : 20 
In an experiment, the two-level system will be influ- 
enced by its environment, which may affect the quantum 
phase of the superposition, alter the effective interaction 
between the levels, or may cause spontaneous decay. For 
qubits in a solid-state environment ) 25 : 26 : 27 : 28 ' 29 all these 
processes may occur simultaneously and hinder qubit ma- 
nipulation. Thus in the context of solid-state quantum 
information processing, a realistic study of qubit manip- 
ulation via Landau-Zener transitions should include the 
influence of environmental degrees of freedom. 

The environment of a quantum system can often be 
described as a bath of harmonic oscillators ] 30 ' 31 : 32 ' 33 : 34 



In some situations, it is known that the dominant envi- 
ronmental effects can best be modelled as a spin bath 
instead ) 35 ' 36 : 37 for example for molecular magnets 3 -^ and 
for Josephson phase qubits 2 — i 28 i 29 at very low tempera- 
tures. 

In the presence of a heat bath, the Landau-Zener dy- 
namics will sensitively depend on the qubit operator to 
which the bath couples^ 8 - Ao and Ramme r 39 : 40 studied 
the LZ problem for the special case in which an ohmic 
heat bath couples to the same operator as the driving 
and derived the transition probabilities in the limit of 
high and of low temperatures. In the limits of very fast 
and very slow sweeps at zero temperature, they found 
that the transition probability is the same as in the ab- 
sence of the heat bath, as was confirmed by numerical 
studies^iil 2 . 

This zero-temperature result was recently proven to 
hold exactly for arbitrary Landau-Zener sweep speeds, 
as a special case of an exact expression for arbitrary 
qubit-bath couplings and spectral densities. 38 Very re- 
cently, Pokrovsky and Sun^ 3 - developed an interesting 
finite-temperature formalism valid for baths with short 
correlation times, in which the exact transition probabil- 
ities of Ref. [38| indeed show up in the zero-temperature 
limit. At sufficiently high temperatures, the qubit expe- 
riences essentially classical Gaussian white noise and the 
Landau-Zener problem can be solved exactly 44 : 45 

A qubit that undergoes a Landau-Zener sweep while 
coupled to another quantum system, is equivalent to a 
multi-level Landau-Zener problem in which two groups 
of levels cross. If all avoided crossings are sufficiently 
well separated, the dynamics consists of effectively inde- 
pendent transitions and one can compute the transition 
probabilities by transfer-matrix techniques^ ' 15 : 46 : 47 In- 
dependent level crossings do not occur for the dissipative 



Landau-Zener problem, for which the adiabatic energy 
spectrum is continuous. Other methods are therefore re- 
quired, for example the exact summation of a perturba- 
tion series^ 

Brundoblcr and Elser— considered a special multi-level 
Landau-Zener problem in which the system starts out 
in the diabatic state whose energy changes faster than 
that of all other diabatic levels. They conjectured that 
the transition probabilities are then given by an expres- 
sion that only contains the velocities of the diabatic lev- 
els and off-diagonal matrix elements of the Hamiltonian. 
The conjecture has long known to be true for the ex- 
actly solvable Demkov-Osherov model, where one level 
crosses several parallel levels . 48 i 49 Analytical proofs of the 
Brundoblcr-Elser conjecture have been given recently for 
several multi-level crossing situations i 50 i 51 i 52 ' 53 i 54 

The physically important situation of a two-level sys- 
tem coupled to a heat bath corresponds to two continuous 
bands of diabatic levels, with the same energy-level veloc- 
ities within each band. A generalization of the theorem 
by Brundoblcr and Elser to this case will be considered. 
This general formula applies to dissipative Landau-Zener 
problems for all those kinds of baths at zero tempera- 
ture for which the initial qubit-plus-bath ground state is 
unique. We discuss the implications of our results for a 
wide range of experiments for which these conditions are 
met. As one of our main results, we present a universal 
aspect of dissipative Landau-Zener transitions, indepen- 
dent of the precise nature of the bath. 

The paper is organized as follows: In Sec. [IT] we review 
the standard LZ formula for a two-level system and gen- 
eralize it to the multi-level case; the explicit calculations 
have been deferred to the Appendix. Section IIIII is de- 
voted to the derivation of the LZ transition probability 
for both a harmonic oscillator bath and a spin bath. In 
Sec. IIV1 we identify a universality in Landau-Zener tun- 
neling that holds true even for baths of nonlinear oscil- 
lators. Finally, we discuss several promising applications 
and experiments in Sec. [V] 

II. LANDAU-ZENER TUNNELING 
PROBABILITIES 

A. Landau-Zener transitions in an isolated qubit 

To set the stage and to introduce our notation, we 
first review the standard Landau-Zener problem for an 
isolated qubit. By "isolated" we mean: not coupled to an 
environment but nevertheless driven by a deterministic 
classical field, whose physical origin will be specified and 
discussed in Sec. [V] The two-level Hamiltonian reads 

/ x y t A ,„ . 

W LZ (t) = -<r z + -a x , (1) 

where <x 2 = | T> <T I - II) Ul and a x = |T)(|| + ||)(T| arc 
Pauli matrices while ||) and |j) denote the so-called dia- 
batic states with the energies ±^vt which cross at t = 0. 



2 



IT) 




time 



FIG. 1: (Color online) Adiabatic (solid) and diabatic (dashed) 
energies for the standard Landau-Zener problem. 

Two parameters determine the dynamics: the constant 
sweep velocity v > by which the energies of the diabatic 
states cross, and the coupling matrix element A between 
these states. Without loss of generality, we assume A 
to be real and non-negative. For A / 0, the diabatic 
states are not eigenstates of the Hamiltonian (Q} , so that 
generally a population transfer is induced. The Hamil- 
tonian is time-dependent and so are its (adiabatic) 
eigenstates. In the limit \t\ ^ A/v, the adiabatic eigen- 
states coincide with the diabatic states. 

The diabatic energies cross, but the adiabatic energies 
±^Vv 2 t 2 + A 2 for A / form an avoided crossing, as 
sketched in Fig. [T] The adiabatic theorem^ states that 
the splitting A prevents transfer of population between 
the adiabatic eigenstates in the adiabatic limit Hv A 2 , 
in other words if the sweep occurs slowly enough. A qubit 
prepared at t — — oo in the initial ground state ||) will 
then end up in the final ground state ||). Beyond the 
adiabatic regime, the dynamics can be rather complex. 
Nevertheless, the population of the diabatic states at t — 
oo can be calculated exactly and is determined by the 
Landau-Zener transition probabilit y 21 i 22 i 23 

P w = | Woo)|i)| 2 = 1 - cxp ( - ^), (2) 

which denotes the probability for a transition to the op- 
posite diabatic state, i.e. the probability to stay in the 
adiabatic eigenstate. Accordingly, = 1 — de- 

notes the probability for a non-adiabatic transition, i.e. 
a jump across the avoided crossing. 

B. Landau-Zener transitions in non-isolated qubits 

We now turn to situations where the qubit is no longer 
isolated. The Hamiltonian of the driven qubit plus its 
environment has the general form 

H(t) = n w {t) + n q -cnv cnv •. (3) 

where Tl onv describes the environment Hamiltonian with 
Hilbert space of dimension M < oo. We assume the 
most general linear coupling between the qubit operators 
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c x , cr y = — i(|T)UI — U)(tl) ; °z and environment opera- 
tors X x,v ' z , in other words we take the qubit-environmcnt 
coupling 

Hq-env = ^ (T V X V . (4) 

^— x,y ,2 

We denote by the eigenstates of the environment 
Hamiltonian H env . 

An important assumption underlying our model §3§ is 
that the qubit-bath coupling (|4|) and the bath itself are 
not affected by the driving. Then at very large times 
t — > ±oo, the qubit Hamiltonian is dominated by the 
time-dependent part, so that all states of the system- 
plus-cnvironment belong to one of two bands: an "up- 
cluster" and a "down-cluster" with energies 

moving upwards and downwards, respectively. 

1. Diabolic basis 

The dissipative Landau-Zener problem is a scattering 
problem in the restricted sense that changes in the qubit 's 
state will occur only during a finite time interval around 
t = 0. In order to exploit this fact that the qubit will 
not flip for sufficiently large \t\, we decompose Tt(t) into 
its diabatic states. These are the eigenstates of the total 
Hamiltonian ([3]) in the the limits t — > ±00. Initially 
and finally, the Hamiltonian is dominated by the term 
proportional to cr z , so that the diabatic basis for the 
qubit is simply given by the states |f) and 

For the environment, by contrast, there is no corre- 
sponding growing energy scale for large \t\. Its diabatic 
states are influenced by the coupling to the qubit and 
depend on the qubit's state. For the "up-cluster" , these 
diabatic eigenstates of the environment are those which 
diagonalize the Hamiltonian projected to the subspace 
|t), i.e. (T|W|T). They are eigenstates of H eav + X z , and 
we denote them by |fc + ) and their energies by £k + . The 
diabatic bath states for the "down-cluster", are de- 
fined likewise so that 

{H cnv ± X z )\k ± ) =e k± \k ± ). (5) 

The diabatic states of the qubit plus the bath read 

|Tfc f ) = |T>|fc+> > (6a) 
U*_) = U)|fc_). (6b) 

where the labels k± = 0, I , 2, . . . are assigned such that 
the energies £& ± are in increasing order, see the sketch in 
Figure [H At asymptotically large times, t — > ±00, the 
diabatic states diagonalize the total Hamiltonian ([3]) and 
hence coincide with adiabatic eigenstates which diago- 
nalize H(t) at a given time t. Note that (J,fc_|tfc+) = 0, 
although in general (fc+|fc'_) ^ Skk' ■ A state of partic- 
ular interest is the adiabatic ground state ||0+) which 
has energy (vt/2 + Sq + ). At zero temperature, it is the 
natural initial state for the Landau-Zener dynamics. 
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FIG. 2: (Color online) Sketch of the diabatic energy levels as 
a function of time for a qubit coupled to a single harmonic 
oscillator. Energies of the states in the "up cluster" increase. 
These states correspond to the qubit state ||). Energies de- 
crease in the "down cluster", where the qubit state is |i). 
According to the "no-go- up theorem" (|AI I [) . the initial state 
|T>0+) evolves to a superposition in which |T,0+) is the only 
"up" state. Energies within a band are separated by the oscil- 
lator energy HQ. For a qubit coupled to an oscillator bath, the 
corresponding crossing clusters would be continuous bands of 
states. 

We now split the Hamiltonian ([3]) into two parts, one 
that is diagonal in the spin index, while the other is off- 
diagonal. The former part consists of all terms propor- 
tional to <r 2 and is diagonal in the diabatic basis ^ . The 
latter part reads 

v = -<r x + ct x x x + cr y xy (7) 

and will be called the bit-flip interaction, since it contains 
all interaction terms of the Hamiltonian §3§ that induce 
a change in the state of the qubit. 

An important feature of the diabatic basis §§§ is that 
all matrix elements of V vanish within each cluster, i.e. 

(Tfc+|v|tfcV) = (ifc_|V|ifcL) = 0. (8) 

This relation will be essential for the application of the 
general formula for the nonadiabatic transition probabil- 
ities as derived in Appendix IA1 

2. Nonadiabatic transition probabilities 

We have now achieved a useful formulation of the dis- 
sipative Landau-Zener problem in terms of two groups 
of diabatic states. If the group of upward moving paral- 
lel levels would consist of merely one state, then tran- 
sition probabilities could be computed with the sim- 
ple independent-crossing formula, for which Brundobler 
and Elser— conjectured that it holds even when succes- 
sive level crossings are not independent. Recent proofs 
show that the independent-crossing formula indeed holds 
exactly, even in more general situationsi 50 i 51 ' 52 i 53 ' 54 As 
stated above, for dissipative Landau-Zener transitions 
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there are two continua of states that cross with constant 
velocity. This physically important situation is addressed 
in the Appendix \A\ where exact nonadiabatic transition 
probabilities are derived in a fairly general setting, with 
the crossing of two continua of parallel states as a special 
case. This setting is a generalization of our recent studies 
in Refs. [Handy! 

For the dissipative Landau- Zener problem, we can de- 
duce the following from Eqs. (jAllj) and (|A21j) : If at 
t = — oo the system starts in a state |T&+) whose dia- 
batic energy is non-degencrate, then the following tran- 
sition probabilities at t = oo are exact: 



P 




27r(Tfc+| V 2 |Tfc+) 

hv 



for k' + = fc^ 
for k', > fca 



For the transition to lower states within the initial band 
of states (k' + < fc + ), we cannot make any statement. 
The second line of Eq. © asserts that states of the 
"up-cluster" that lie above the initial state will finally 
be u npop ulated. This no-go theorem was formulated in 
Refs. l56l and l57l and we think that it is more aptly de- 
scribed by the name "no-go-up theorem" . 

A case of particular interest is that of the initial state 
|T0+), which is the ground state of the entire system. 
For all bath models employed below, the ground state 
is unique, so that relation ^ applies. Then final states 
with k' + < k + do not exist, while the occupation of final 
states with k' + > k + is forbidden by the no-go-up theo- 
rem. Thus, provided that the final qubit state is ||), the 
environment will end up in its ground state. 

It is the final transition probabilities P|->-f and Pf— >.j, 
for the qubit that interests us most, irrespective of the 
final state of the environment. By tracing out the envi- 
ronment, i.e. by performing the sum over k' + , we find 



T->T = ex P 



ttW 2 
2hv 



= l-P 



T->1 



with the ground-state expectation value 
W 2 = 4(t0+| V 2 |T0+) 



(10) 



(11) 



These are the two central equations for dissipative 
Landau-Zener transitions at zero temperature. The 
ground-state expectation value W 2 formally replaces 
squared the tunnel matrix element A 2 in the original 
Landau-Zener expression ([2]). 



III. DISSIPATIVE LANDAU-ZENER 
TRANSITIONS IN VARIOUS ENVIRONMENTS 

We are now in the position to study Landau-Zener 
transitions for specific baths with linear couplings to the 
qubit. We will focus on the two most important model 
baths in quantum dissipation, namely a bath of harmonic 
oscillators and a spin bath. In both cases, we will restrict 
ourselves to zero temperature, so that the natural initial 



state is the diabatic ground state |T0+) of the system plus 
the environment. Both for the harmonic-oscillator bath 
and for the spin bath, the ground state is non-degenerate, 
so that formula (JTDJ) can be applied. The essential steps 
that remain are first to identify and characterize the dia- 
batic ground state |T0+)> and then to compute the expec- 
tation value Applications to specific experiments 
will be discussed in Sec. IVl 



A. Harmonic-oscillator bath 

We first consider the case in which a qubit interacts 
with a standard bosonic bath consisting of harmonic os- 
cillators. The Hamiltonian is as in Eq. Q, with the 
environment Hamiltonian 



N 



(12) 



consisting of N harmonic oscillators with frequencies flj . 
Zero-point energies do not play a role here and are ig- 
nored. The &j and bj denote the usual creation and anni- 
hilation operators of the oscillator j. We leave the num- 
ber of oscillators N finite at first but eventually take it 
to infinity in a continuum limit. Furthermore, we assume 
the most general linear qubit-oscillator coupling 



N 



where the second sum specifies the environment oper- 
ators X v defined in Eq. ((3]). Since the coupling (fT5|) 
also includes the ct v interaction, this constitutes a gen- 
eralization of the spin-boson model that we considered 
in Ref. ,38,. The parameters jj determine the coupling 
strengths, while the parameters X" define the "coupling 
directions" and are conveniently expressed by the spher- 
ical coordinates J?, and a?, as 



X^ = sin i3j cos ipj , 
Aj = sin $j sin ipj , 
Xj = cos'dj. 

The bit-flip interaction (O then becomes 



(14a) 
(14b) 
(14c) 



A N 

1 3=1 



7jSini?j 



(cr x cos (fj + <j y sin y>j ) (bj + ) . 

The diabatic states of the environment are determined 
by diagonalizing the Hamiltonian 

n cnv + X z = E Mljb]bj + E ^ cos ^ (&} + bj) (16) 



E 



(17) 
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In order to obtain the diagonal form in the last line, we 
introduced the shifted annihilation operators 



J 3+ 



7jCQs^ 

2m -i 



(18) 



and the reorganization energy of the jth oscillator, Ej = 
(7 2 /4M2j) cos 2 The latter quantity denotes the en- 
ergy shift of the oscillator ground state owing to the cou- 
pling to the qubit. From Eq. (|17p. it becomes immedi- 
ately clear that the diabatic ground state has to fulfill 
the relation 



o. 



(19) 



while the excited diabatic states are created by applying 
the operators to this diabatic ground state. 

We now write the bit-flip interaction V with the shifted 
operators bj, and employ relation (| 19|) to evaluate the 
ground state expectation value (fTTj) . After some algebra, 
we obtain^ 



W 2 = 



3 



7 sin <d 



J ! 



(20) 



which allows one to compute the Landau-Zener transition 
probabilities (fTU|) . Note the (^-dependence, generalizing 
our recent work^ which shows that for A ^ 0, it makes a 
difference for Landau-Zener transition probabilities what 
types of off-diagonal qubit-oscillator couplings exist. 



1. Identical coupling angles 

For system-bath models, it is frequently assumed that 
all bath oscillators couple to the central quantum system 
via the same coordinate. In our model, this corresponds 
to the case in which all coupling angles are identical, 



0, 



d and <pj = <p. Then Eq. ([20]) becomes 



W\$,<p) 



-£ sin(2i?)e- l¥ T + Ssm 2 d, (21) 



in terms of the integrated spectral density^ 



5' 



v 2 = i- f 

4 Wo 



and the energy 



-in 



duo J(lo) 



(22) 



(23) 



which equals four times the total reorganization energy 
E = Y^,j Ej . These quantities S and Eq were presented 
as frequency integrals over the spectral density J(u>) of 
the bath, with the latter defined as 



N o 2 



(24) 
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FIG. 3: (Color online) Landau-Zener dynamics for a qubit 
with A = 0, in all three cases shown off-diagonally coupled 
via (j x to three oscillators. The various oscillator frequencies 
Qj are given in units of s/v/H. All coupling strengths have the 
same value jj — y/hv/3. The dotted line marks the analytical 
final transition probability corresponding to Eq. Ij28fl . 



In a continuum limit, the spectral density becomes 
a smooth function of frequency. At low frequencies, 
one typically observes a power-law behavior for J(w), 
whereas the qubit becomes insensitive to high frequen- 
cies, as characterized by a cutoff frequency lj c . An im- 
portant class of spectral densities is therefor e) 31 ! 34 



S —-Uj/uJr 

— auj e 1 



(25) 



For such spectral densities, one immediately obtains the 
two global quantities whereby the bath influences LZ 
transition probabilities: 



„ ah „„. , 
E Q = ^ C T( S ), 



(26) 
(27) 



where r(x) denotes the Euler Gamma function. 

a. Off- diagonal coupling. For •& = tt/2, the qubit 
interacts via its off-diagonal operators er x and <j v with 
the environment, whereas the LZ driving affects the qubit 
only via <j z . Equation (|2ip becomes 



W 2 (Tr/2,tp) = A 2 + S. 



(28) 



Interestingly enough, the Landau-Zener tunneling proba- 
bility is then fully determined by the integrated spectral 
density S. In particular, there is no dependence on the 
oscillator frequencies Qj . This result is nicely illustrated 
in the simple example of Figurc[31 showing Landau-Zener 
dynamics of a qubit that is coupled to only three oscil- 
lators. The oscillator frequencies are varied, while the 
qubit-oscillator couplings are kept constant. The dy- 
namics at intermediate times depends on the oscillator 
frequencies, but the final transition probability does not. 
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Note that FF 2 (7r/2, tp) in Eq. (gSJ is independent of tp, 
so that in case of only off-diagonal coupling the relative 
weight of the interactions via ar x and er y drops out of 
the final occupation probability. Since S > 0, it is clear 
from Eqs. (fTTj) and (|28|) that a off-diagonal coupling al- 
ways enhances the occupation of the final ground state 
||) as compared to the case without dissipation. This has 
an intuitive explanation: the zero-temperature oscillator 
bath partially succeeds in cooling the qubit down to its 
instantaneous ground state at any time during the level 
crossing. 

b. Diagonal coupling. For — 0, the qubit interacts 
with the environment through the diagonal Pauli matrix 
<j z . This interaction induces pure dephasing between 
the states \\) and |J.). This driven spin-boson model 
with diagonal interaction has so far been the standard 
model for discussing Landau-Zener transitions in dissipa- 
tive environmentsj 39 ' 41 i 44 ' 45 i 59 For diagonal coupling, the 
bit-flip interaction ([15]) simply becomes V = (A/2)a x , so 
that 
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FIG. 4: (Color online) Transition probability P-f— ►! as a func- 
tion of the internal coupling A for Eq = 2\fhv and S = O.bhv. 
The angles i? = and i? = 7r/2 correspond to purely diagonal 
and off-diagonal coupling, respectively, for which the proba- 
bility is independent of tp. 



W(0,p) 2 = A 2 , (29) 

and the Landau-Zener transition probability (|I0[) coin- 
cides with the standard expression Eq. ^ for an isolated 
qubit. This bath independence of the transition probabil- 
ity for diagonal coupling was predicted by Ao and Ram- 
mer in the limits of fast and slow Landau-Zener sweeps^ 9 - 
Here we find that at zero temperature, it holds exactly 
for all diagonally coupled harmonic-oscillator baths, and 
for all coupling strengths A and sweep velocities vM- 

c. General coupling. When the oscillators neither 
couple purely off-diagonally (t? = tt/2) nor purely diag- 
onally (d = 0), the Landau-Zener probability generally 
exhibits a non-monotonic dependence on the tunnel cou- 
pling A. This is shown in Figure U for various angles 
p and $. Most interesting is the comparison to the non- 
dissipative case, which as we saw coincides with the result 
for diagonal coupling ($ — 0): Any dissipative Landau- 
Zener probability lower than the curve for d = is coun- 
terintuitive. Such situations occur, however: for several 
values of $ and for a sufficiently large tunnel splitting 
A, the bath coupling reduces Pf_>j(oo), i.e. dissipation 
enhances the population of the final excited qubit stated 

This counterintuitive behavior is most significant for 
the angles <p — 0, = tt / '4 and when the squared re- 
organization energy is large compared to the integrated 
spectral density Eq ^> S. This means that the coun- 
terintuitive behavior would be observable only in rather 
exceptional situations, which did not shape our intuition. 
For example, for the spectral density (|25|) . a significant 
reduction of the ground state population by increase of A 
requires a very strong qubit-bath coupling a > 1. In the 
opposite limit a -C 1, which is relevant for quantum in- 
formation processing, the bath for all practical purposes 
induces the more intuitive tendency towards the ground 
state as A is increased, as we found above for purely 
off-diagonal or purely diagonal couplings. 



B. Spin bath 

Let us now turn to the case in which a qubit interacts 
with an ensemble of otherwise non-interacting two-level 
systems forming a spin bathi 35 i 36 i 37 The total Hamil- 
tonian is again assumed to be of the general form (|3j). 
with the standard Landau-Zener Hamiltonian ([1]) for the 
qubit, but now with the bath Hamiltonian 

N 

3=1 v=x,y,z 

where are the Pauli matrices for the jth bath spin. 
The most general linear qubit-bath coupling reads 

N 

H q _env = E ^EtW' ( 31 ) 
iy—x.y^z 3 — 1 

where the second sum defines the operators X v as a lin- 
ear combination of the spin operators r^, with coupling 
constants 7J. The bit-flip interaction V then becomes 

A N 

V 7 <r, • X (tJ^tJ + ^ayii) . (32) 
3=1 

As for the bosonic bath in Sec. IIII A[ we wish to apply 
the general formula (fl"C)|) for the transition probability. 
For that, we need to determine the diabatic eigenstates 
of the qubit plus the spin bath. For large \t\, the time- 
dependent term vta z /2 dominates and, therefore, pro- 
vides the diabatic qubit states |f) and ||). Consequently, 
the diabatic spin-bath states are determined by the op- 
erator 

(s\H(t)\s) =J2Kin,p s = t,4, (33) 
3 
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where "+" refers to s = f and "— " to s = J., while 

n% inj = ±^T Z + J2 (34) 



v— x,y .z 



determines the state of the jth bath spin. The Hamil- 
tonian (|34|) is readily diagonalized and its eigenenergies 
— Aj -t and Aj t ± are determined by 



(B* ±7-) 2 - 



(35) 



For the evaluation of the Landau-Zener transition prob- 
ability we will not need an explicit expression for the 
ground states \0j,±) of Hf pin j. It suffices to know 
that the ground states satisfy the eigenvalue equation 



n fpinj l°.?±) = -Aj± \0j±) . Consequently, 



(0j+W|0 i+ ) = -Bf/A j+ , (36a) 
(0 J+ kJ|0 J+ ) = -B*/A j+ , (36b) 
(0 J -+W|0 i+ ) = -(B*+ 7 J)/A i+ . (36c) 

Since the bath spins do not interact with each other, the 
diabatic ground state |0±) is the direct product of the 
states |0±). 

With this ground state defined, we are now in the po- 
sition to employ formula (|10|) . Inserting relations 
into (fTTj) . we obtain 

IF 2 = A 2 -4A£( 7 *) 2 ^ + 4 ]T[( 7 J) 2 + ( 7 J) 2 ] 



B v -B v ., 



'A 7 ' 7 ' A i+ A iH 7 ^ Aj+Aji 



(37) 



which determines the Landau-Zener transition proba- 
bility (|T0|) . The last term stems from the commuta- 
tor [cr x ,(Ty\ = 2icr z . Recall that the general Landau- 
Zener formula (|10p was derived under the assumption 
that the diabatic qubit-plus-bath ground state at t = — oo 
is non-degenerate. Therefore, our results do not ap- 
ply to parameter sets for which any bath spin obeys 

B ! = B j = ( B i + % 2 ) = so that A J+ = °- 

a. Diagonal basis for spins. First, we consider the 
physically important special case that all and -BJ van- 
ish, so that the eigenvalues become A J+ = \B? + 7j|- The 
ground-state expectation value W 2 which determines the 
tunnel rate then assumes the more compact form 



n-^ = A-+4f;(^ + .s- / r;'] 

3=1 



(38) 



where Sj = sign(_B| +7*)- We can see that in this 
special case, the transition probability is independent of 



the eigenfrequencies Aj/K; the probability essentially de- 
pends on the qubit-spin coupling strengths 7 J. This fre- 
quency independence resembles the transition probabil- 
ity l|28p for the off-diagonally coupled bosonic oscillator 
bath, although the qubit-spin coupling in the present case 
is not necessarily off-diagonal. 

b. Pure dephasing. Interestingly, it follows from 
Eq. (J38]) that if all (7J + stf?) vanish, that then the 
tunneling probability P[—t± equals the standard LZ prob- 
ability despite coupling to the bath. This includes 
the important case of pure diagonal coupling for which 
7J = 7J = 0. This bath independence confirms and 
generalizes the very recent perturbative calculations by 
Wan et al£2. for dissipative LZ transitions in a spin bath, 
where it was assumed that A is small. We find the same 
exact bath-independent transition probability in case of 
a diagonally coupled spin bath as we found earlier in 
Sec. IIII Al for a diagonally coupled bosonic bath. This 
striking result is discussed further in Sec. IIV1 

c. Robustness under dephasing. Without loss of 
generality, the energies B* of the spins can be chosen 
positive. For spin baths away from the very-strong cou- 
pling regime, the coupling constants 7J will be much 
smaller than the corresponding energies B?, so that 
sign(-B| +7|) = 1. Thus unless the qubit-spin coupling- 
is very strong, the LZ transition probability correspond- 
ing to Eq. (|3"5|) is independent of the 7J. Interestingly, 
we find that this holds true even when 3> 7j :!y , i.e. if 
dephasing is much stronger than relaxation. 

Moreover, in many models for decoherence by spin 
baths, the coupling of the bath to the cr v operator simply 
does not occur, in which case Eq. (j3"5|) reduces to 



W 



A 



N 



3 = 1 



(7i 



(39) 



Note that this result is obtained without making any as- 
sumptions on the dephasing strengths 7?. In other words, 
the LZ transition probability only depends on the inte- 
grated spectral density for relaxation, and is fully inde- 
pendent of the dephasing strength of the spin bath. As 
discussed in Sec. [Vl this is important in experiments. 



IV. UNIVERSALITY OF BATH-INDEPENDENT 
NONADIABATIC TUNNELING PROBABILITY 



As was found in the previous section, when the qubit 
interacts with a bosonic bath or with a spin-bath only via 
the operator cr z , and the total system starts in its ground 
state, then the probabilities for dissipative Landau-Zener 
transitions coincide with the standard tunneling proba- 
bility ([2]). Thus, a natural question arises: Is this a co- 
incidence, a speciality of the two baths that we studied, 
or does it hold more generally? We will now show that 
this holds in general, regardless of the specifics of the 
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environment. Consider the Hamiltonian 

H = y(T 2 + yITj; + (T Z X Z + H c 



(40) 



which is the general Hamiltonian Q specified for diago- 
nal coupling (X x = %y = 0). As before in Sec. ITTBl the 
operator X z is the environment operator with which the 
environment is coupled to the qubit, and 7i nv is the en- 
vironment Hamiltonian. Further specifications of these 
operators need not be given for our reasoning. Note that 
the qubit-bath interaction er z X z does not commute with 
the standard LZ Hamiltonian at any time if A ^ 0. Tak- 
ing this at face value, it is tempting to assume that the 
transition probability P-[—>i will be affected by a bath 
that causes pure dephasing. Our analysis, however, will 
reveal that this is not the case. 

As before, a complete set of diabatic qubit-plus-bath 
states \1k + ) and can be found, where the shifted 

diabatic bath states |fc±) are eigenstates of (Tt en v ±X Z ). 
We assume that (H. cm + X z ) has a nondegenerate ground 
state |0+). Now since the bath does not induce bit flips, 
the bit-flip operator ([7]) simply reduces to the internal 
bit-flip interaction, i.e. V = Acr x /2. Hence W 2 = A 2 . 
Thus from Eq. (fit?)) we find the standard Landau-Zcncr 
transition probability for a qubit diagonally coupled to 
an arbitrary bath at zero temperature, 



f-r-a = 1 - exp - 



7rA 2 
2hv 



(41) 



It is truly remarkable that diagonal coupling to the en- 
vironment does not affect the final transition probabil- 
ity, whatever the nature of the environment and however 
strong the coupling operator X z may be. In other words, 
this zero-temperature transition probability for diagonal 
bath coupling (|4Tj) indeed holds universally. It may have 
been simple to derive it from Eq. (|10|). but the physical 
implications are most important: Landau-Zener transi- 
tions are insensitive to pure dephasing at zero temper- 
ature, irrespective of the nature of the bath or of the 
bath-coupling operator X z . 

To illustrate the universality that we just found, let us 
now also consider the coupling to a collection of nonlinear 
oscillators 



N 

i=i 



Pi 
4! 



(42) 



that couples to the qubit via the interaction operator 



N 

i=i 



(b] + b 3 ). 



(43) 



LZ sweeps pas resonances of nonlinear oscillators are of 
practical interest, since nonlinear oscillators are currently 
used for the readout of flux qubits?^ For a numerical test 
of the predicted final transition probability, we take the 
situation in which the qubit is diagonally coupled to two 
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FIG. 5: (Color online) Time evolution of spin-flip probability 
for a qubit with A = 0.5\/fi/u diagonally coupled to two har- 
monic oscillators, two nonlinear oscillators, and seven spins, 
respectively. The harmonic oscillators are specified by Qi = 
Q.ly/v/h, Q2 = 0.5-\/vJh,, 71 = 2\fEv, and 72 = Qy/hv, while 
the nonlinear oscillators, in addition, have j3\ = fl% = 
The values of the B" and the 7J are randomly chosen from the 
range [— VHv/10, -\/fi/u/10]. In all three cases, the transition 
probability converges to the universal value (14111 . 



of these nonlinear oscillators. Figure [5] shows the cor- 
responding time-evolution of the probability P-f_»^(t) for 
the qubit to be in the "down" state. It also shows the 
dynamics in case the qubit couples to two linear oscilla- 
tors with the same parameters, except that now f3\ t 2 — 0. 
Furthermore, the effect of a diagonally coupled spin bath, 
a special case of Eq. (l30|) . is also shown. We consider the 
case in which the spin bath consists of seven spins. The 
internal interaction A has the same value for all curves in 
Fig. [5] As shown in the figure, at intermediate times the 
transition probability shows variations depending on the 
specifics of the environments, but the final probabilities 
indeed all converge to the universal value (HTj) . 



V. EXPERIMENTAL RELEVANCE 

Landau-Zener transitions are commonly used in exper- 
iments to determine the interaction A between diabatic 
energy levels^ One usually takes the standard LZ for- 
mula (|2|) as a starting point. For fixed sweep speed v of 
the energy levels, the only unknown in the equation is 
the internal interaction A. By varying v, one can deter- 
mine the Landau-Zener probability -Pf— >f and, hence, the 
A. This in turn is a method to test the validity of the 
two- level result ©. 

However, if the qubit is coupled to an environment that 
causes relaxation, then Eq. (|28p shows that one actu- 
ally measures VA 2 + S for a qubit in an oscillator bath, 
where S denotes the integrated spectral density. Vari- 
ation of the sweep speed v does not help in extracting 
A and S separately from a Landau-Zener experiment. If 
the bath only couples diagonally and thereby causes only 
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pure dephasing, then at zero temperature, the bath does 
not influence the Landau-Zener tunneling probability at 
all. 

In most situations, the environment will cause both re- 
laxation and pure dephasing. We find that dephasing will 
hardly change the relaxation-dependent transition proba- 
bilities, see Eqs. (|2ip and (|38|) . In many experiments pure 
dephasing is much faster than relaxation, in other words 
T 2 * times are much shorter than T\ times. It will therefore 
often occur that S sin 2 $ <C A 2 and, hence, relaxation can 
be neglected on the time scale of the Landau-Zener tran- 
sition. Then our results imply that with a Landau-Zener 
experiment, one can accurately determine the internal 
interaction A even in the presence of strong dephasing, 
such that iScos 2 d 3> A 2 . Landau-Zener experiments are 
therefore a surprisingly reliable tool to determine tunnel 
splittings. We now consider some specific applications. 



A. Circuit QED 

Circuit QED 62 i££ is a superconducting-circuit analogue 
of optical cavity QED: A charge qubit^ 2 - or a flux qubit 63 
is coupled so strongly to a a quantized harmonic oscilla- 
tor in the circuit that Rabi oscillations can be observed in 
a solid state environment. Because the circuit-QED sys- 
tem is so highly tunable, it enables the study of quantum 
dynamics of open quantum systems in new parameter 
regimes. Recently, we proposed to use LZ transitions in 
circuit QED to generate single microwave photons^ by 
choosing parameters such that the qubit-oscillator cou- 
pling is off-diagonal and the internal interaction A van- 
ishes. Here we consider an arbitrary A and discuss the 
possibility of qubit-oscillator couplings other than off- 
diagonal. 

The Hamiltonian describing a charge qubit interacting 
with the transmission-line resonator is 



H(t)=H q (t)+H q . osc (t)+H Q 



(44) 



where the different terms describe the qubit, the oscil- 
lator, and the qubit-oscillator coupling, respectively. If 
the dynamics is essentially restricted to two states with 
N and N + 1 Cooper pairs in the box, then the Hamil- 
tonian becomes beautifully simple^ and in our notation 
reads 



H a — — 



Hn 



2 



A 

2 



7 2 (l-27V K dc ) 2 



H 



7(& + 6 t )<r x , 



(45a) 

(45b) 
(45c) 



with the charging energy Eq, the oscillator frequency u T , 
and the coupling strength 7. We presented the Hamil- 
tonian in so-called tunneling representation, which is the 
basis in which the charge states \N) and \N + 1) are 
eigenstates of <r x . Since the Josephson link is imple- 
mented by a dc SQUID, the Josephson energy Ec(t) = 



-Ej,max cos[$(i)/ ( l > o] can be manipulated upon variation 
of the flux that penetrates the SQUID^ where 4> 
denotes the flux quantum. The dimensionless quantity 
Ng c is proportional to an applied gate voltage which can 
be used to adjust the internal coupling, which in terms of 
the control parameters reads A = 4(1 — 2N^ c )(j 2 /huj,- — 
E c ). 

Since the internal coupling A is tunable for this setup, 
one can use Landau-Zener transitions to determine the 
integrated spectral density S for the purely off-diagonal 
coupling (|45c[) . In the one-oscillator model (|4"5| , S — 
7 2 , so that W 2 = A 2 + 7 2 . Thus, by Landau-Zener 
sweeps past the oscillator resonance for different values 
of A, Eq. (|2"8|) allows one to determine the qubit-oscillator 
coupling strength 7 2 in independent ways. 

The Hamiltonian (145[) features a possibly nonzero in- 
ternal interaction A between the qubit levels, which 
was assumed zero in the setting of this superconducting 
circuit for which single-photon generation was recently 
proposed^ Besides the transitions |f0) — > |J,1) that gen- 
erate a single photon, for A ^ there is now also the 
process |f0) — > 1 4,0) which flips the qubit without chang- 
ing the cavity state. Therefore, a reliable single-photon 
generation requires that the transitions induced by A are 
not relevant, while the qubit-oscillator coupling has to be 
sufficiently large?^ such that 7 2 3> Hv. Fortunately, as 
long as A < Hu> T , efficient single-photon generation is 
very well possible by sweeping the qubit energy only on 
a finite frequency interval small compared to w T but large 
compared to the width of the cavity resonance. 

It would be interesting to test in circuit QED our 
predictions for other than off-diagonal qubit-bath cou- 
plings. This requires magnetic fields and gate voltages 
to be swept simultaneously. The problem with sweeping 
gate voltages, however, is that this displaces the ground 
state of the oscillator,— which practically means that A 
cannot be used as a time-dependent control parameter. 
Consequently, the electric-plus-magnetic sweeping in cir- 
cuit QED falls outside the class of models that we con- 
sidered. It seems that Landau-Zener transitions with a 
general qubit-oscillator coupling angle d as in Eq. ([±3]) 
cannot be engineered in state-of-the-art circuit QED. 

We already mentioned the width of the cavity reso- 
nance, but actually in the model (|4"S")) , we assumed that 
the quality factor of the resonator was infinite. The spec- 
tral density of the oscillator bath then consists of a sin- 
gle delta peak at the cavity resonance frequency, so that 
S = 7 2 . For the non- ideal resonance, the spectral density 
become a 25 i 26 



J{u>) = 



(UJ — UJ r ) 2 + K 2 LU 2 (V 2 



(46) 



where the dimensionless parameters a and k measure the 
qubit-oscillator strength and the width of the resonance 
peak, respectively. For small frequencies w <C uj r , the 
spectral density J(lu) ~ anu>, as for an ohmic bath. Now 
suppose that we let the qubit undergo a LZ sweep from 
frequency zero across the resonator frequency until a fre- 
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quency ui max > (1 + k)uj t . Now since for frequencies 
\uj — uj r \ > Kuj r the spectral density falls off rapidly, a hy- 
pothetical continuation of the LZ sweep from frequency 
w max to infinity would hardly change the transition prob- 
ability. Therefore, by combining Eqs. ([22]) and ([46]) we 
find the integrated spectral density 

= 0(^1 j / 2-£ \ «\ ofrrf 

4^ \ \kVA^J 2] 4 

(47) 

Inserting this expression into Eq. [)28p gives an accurate 
value for the cavity-induced LZ transition probability. 
The last identity in (|47ll holds in the experimentally re- 
alized good-cavity limit k«1. In this limit one can ap- 
proximate the spectral density (J3H1) by a single oscillator 
with qubit-oscillator coupling 7 2 = a(liw r ) 2 /4, reproduc- 
ing the model (|4"B")l . 

This leads us to the important conclusion that LZ tran- 
sition probabilities can be computed exactly for qubits 
swept past narrow or broad cavity resonances alike. The 
total strength S of the atom-cavity coupling determines 
the LZ transition probability, and for weak dissipation S 
is independent of the scaled cavity width n. Dissipative 
Landau-Zener transitions are also robust in this respect. 



B. Cavity QED and photonic crystals 

Atoms couple off-diagonally to the electromagnetic 
field, which for many practical purposes can be consid- 
ered as a bath at zero temperature. Atomic energies 
can be swept by applying dc electric and magnetic fields, 
which give rise to Stark and Zeeman shifts, respectively. 
Atomic resonances are usually narrow but hard to sweep 
by large percentages. Artificial atoms are more tunable, 
but their resonances are also broader. Our results ap- 
ply to situations where resonances can be swept by much 
more than their width. 

An atom in free space feels a spectral density that is 
quadratic in the frequency, J(u>) cx (uj/ujq) 2 , with no cut- 
off frequencies in the optical regime, so that S — > oo. 
During a Landau-Zener sweep of an atom in free space, 
the atom will finally have decayed spontaneously to its 
ground state. However, the spectral density felt by the 
atom, also known as its local optical density of states^ 6 - 
can be engineered by changing its dielectric environ- 
ment. In photonic crystals for example, the periodicity 
of the refractive index on the scale of an optical wave- 
length Ao = 2nc/u)o may create a photonic band gap of 
width. Within this gap, the spectral density J{lo) ide- 
ally vanishes and in practice it can be strongly reduced, 
whereby spontaneous emission by the atom is strongly 
suppressed^ 

By making a controlled point defect or line defect de- 
fect in the vicinity of the atom that breaks the peri- 
odicity of the photonic crystal, a narrow defect mode 
may be created^ 8 - within the spectral gap, as sketched 




FIG. 6: (Color online) Sketch of spectral density for an atom 
near a local defect in a photonic crystal with a band gap. The 
quadratic free-space spectral density is modified by the crystal 
that creates a spectral gap around ujq. A narrow defect mode 
inside a broader band gap allows a controlled atom-defect 
interaction via LZ sweeps of the atomic transition frequency 

LUA(t). 

in Fig. [HI Ideally, this would allow cavity QED exper- 
iments to be performed within a photonic crystal, and 
progress is made in this direction] 67 ! 68 We propose to do 
LZ sweeps of the atomic frequency u>a(£) around the de- 
fect frequency but within the band gap. This will allow 
the creation of atom-defect entanglement and of single 
photons in the defect mode, in quite the same way as in 
circuit QED. 18 

C. Molecular nanomagnets 

The energy levels of molecular nanomagnets can be 
swept by switching on dc magnetic fields. 3 Higher-excited 
states have higher magnetization and these are excited 
more when switching rates are high. These nanomag- 
nets are not intrinsic two-level systems and indeed many 
successive LZ transitions are observed upon increasing 
the magnetic field. They can be cooled until tunnel- 
ing rates become temperature independent. 3 In recent 
experiments^ this happens rather abruptly at a tempera- 
ture of 0.6 K. In this low-temperature pure quantum tun- 
neling regime, our predictions apply. LZ transitions are 
commonly used in experiments to determine level inter- 
actions between energy levels of molecular nanomagnets. 
Energy relaxation due to thermal environment usually 
changes the effective energy gap. 69 But as we showed, LZ 
transitions for qubits are robust under dephasing, which 
on a qualitative level agrees with previous lab experience 
about LZ transitions in nanomagnets 

Molecular nanomagnets show many avoided crossings 
as the magnetic field is varied. These crossings are often 
well separated, so that one can ask what type of bath in- 
fluences the effective-two-qubit dynamics around a par- 
ticular level crossing. Neither the bath nor its coupling 
to the qubit are known precisely, but at very low temper- 



11 



atures the main source of decoherence will stem from the 
coupling to other electronic and nuclear spins, so that 
a spin bath seems the appropriate model. Our results 
of Sec. IIII Bl show in detail how transition probabilities 
depend on both the internal interaction A and the in- 
tegrated spectral density S. Although strong dephasing 
is probable in nanomagnets, our analytical results prove 
that this hardly affects the LZ transition probabilities 
in the experimentally accessible regime where tunneling 
rates are temperature independent A 

A theory for multiple dissipative LZ transitions in 
molecular magnets was developed by Leuenberger and 
Loss.' This theory presumes that in between successive 
LZ transitions, all quantum coherences have been washed 
out due to dephasing while dephasing during the indi- 
vidual LZ transitions does not influence the transition 
probabilities. This assumption is rather natural if the 
LZ tunneling time tlz is much shorter than the pure de- 
phasing time, but obviously requires proof for slower LZ 
transitions. This is an important application of our re- 
sults: our derivation of the universal bath independence 
of individual LZ transitions under dephasing proves that 
the theory of Ref. UM is more widely applicable than one 
would have guessed previously. 

D. Quantum computer as a spin-bath simulator 

An ideal quantum computer is a collection of qubits 
whose energies and internal interactions are tunable, and 
whose mutual interactions are tunable as well. One 
method for single-qubit manipulation is a Landau-Zcncr 
sweep. This has been realized in recent experiments on 
superconducting qubits i 13 ' 17 ! 20 If the interactions with 
the other qubits are not exactly zero, then our spin-bath 
result ([57)1 predicts an effect on the "single-qubit" LZ 
transition probability. This way, the LZ sweep provides 
a test of the settings of the quantum computer in the 
operational space of the qubit. 

On the other hand, our predictions can be tested very 
carefully by controlled variation of qubit-qubit interac- 
tions in the quantum computer. Indeed, the quantum 
computer could be seen as a quantum simulator, in our 
case of the effects of a spin bath on the LZ transi- 
tion probability. A full-fledged quantum simulator does 
not exist yet, but recent experiments on a system of 
four superconducting qubits with tunable couplings look 
promising^! 

E. Superconducting qubits and their spin baths 

Josephson phase qubits at very low temperatures 
(20 mK) exhibit decoherence mainly due to interactions 
with two-level microwave resonators.— Since it is found 
that qubit losses strongly depend on driving amplitudes, 
these resonators cannot be described as a bosonic bath. 
Rather, they are nowadays thought to be charge two-level 



systems. 29 The existence of these resonators has impor- 
tant implications for the operation of superconducting 
qubits. Spectroscopic measurements have shown that at 
fixed energy, the qubit often resonantly interacts with 
only a single resonator. Moreover, coherent quantum os- 
cillations between a qubit and such a single resonator 
were observed. 2 - Remarkably, decoherence times of these 
microscopic two-level systems are larger than that of the 
qubit. This in turn has led to the very recent proposal 
to use instead these microscopic resonators as qubits for 
quantum information processing. 72 

Our exact result (f3"7) for the Landau-Zener transi- 
tion probability of a qubit in a spin bath in the low- 
temperature tunneling regime can be an important tool 
for analyzing further the properties of the microscopic 
resonators and their couplings to the qubit. Our assump- 
tion of a spin-star configuration, i.e. a bath of mutually 
non-interacting spins, is probably correct for the spin 
bath of the phase qubits 3^ Monte-Carlo simulations in- 
dicate that narrower qubit-spin resonances are shadowed 
by larger onesi^ 9 - For that reason it is important that for- 
mula (|37p holds generally, whether the qubit resonantly 
interacts with one microwave resonator at a time or not. 

During a so-called fast-pulse measurement of the state 
of the qubit, the qubit energy moves in and out of res- 
onance with many of these resonators^ As a conse- 
quence, the resonators reduce the fidelity of the measure- 
ment. Actually, this effect of the fluctuators has already 
been estimated in terms of multiple LZ transitions in 
Ref. There it was assumed that the resonators couple 
off-diagonally via Ylj=i lj CT x T i> with, coupling strengths 
ilj ) 2 given by the size of the splittings as measured in the 
qubit spectroscopy. Our results (|37D and especially ([55]) 
show that the formula used in Ref. [28| becomes exact for 
a zero-temperature spin bath that couples off-diagonally, 
even in case of overlapping spurious resonances ^ It is 
unfortunate that microscopic resonators reduce the fi- 
delity of the fast-pulse readout method, but it is good 
to know how much. More good news is that the same 
fidelity reduction would be obtained even if there would 
be additional dephasing by the spins. 

Landau-Zener sweeps for phase qubits can be relevant 
for one more reason)^ the precise spectral distribution 
of the spins will be sample dependent. Moreover, the 
distribution for a single sample varies on a time scale of 
daysi 7 - 2 - Fortunately, the LZ sweep measures a "global" 
property of the spin bath, namely its integrated spectral 
density S. It is fair to assume that S will vary less from 
sample to sample and from one day to the next than 
J(ui) at a fixed frequency. We therefore suggest to use 
LZ sweeps as a robust way of "gauging" and characteriz- 
ing the spin-bath environment of superconducting phase 
qubits. 
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VI. SUMMARY AND CONCLUSIONS 



We studied the effect of various zero-temperature envi- 
ronments on the Landau-Zener transition probability of a 
qubit. The main result of this paper is the corresponding 
generalization of the bit-flip transition. The mathemat- 
ical form of this result, just like the standard Landau- 
Zener formula, is charmingly simple: from the total 
Hamiltonian one identifies the bit-flip operator V of the 
qubit and computes the expectation value W 2 = 4(V 2 ) 
for the initial diabatic ground state of the qubit plus its 
environment. The transition probability follows formally 
by replacing in the standard expression for the two-level 
LZ problem the factor A 2 by W 2 . Following this recipe, 
we have calculated in Sec. Mil the LZ tunneling probabil- 
ities for oscillator baths and for spin baths, which repre- 
sent the two most important environments in quantum 
dissipation research. These examples illustrated that dis- 
sipative Landau-Zener transition probabilities in general 
depend on the type of environment and on the way the 
qubit couples to it. 

However, when the qubit-bath coupling is of the di- 
agonal type, causing so-called pure dephasing, then the 
tunneling probability coincides with the original LZ prob- 
ability, regardless of the details of the environment. We 
expect that this universal behavior is observable and 
important for a wide variety of applications. A bath- 
independence of another kind was found for a qubit swept 
past a broadened (circuit) cavity resonance: the tran- 
sition probability turned out to be independent of the 
quality factor of the cavity. Since transition probabili- 
ties are also independent of the environment parameters 
in some phenomenological models with non-Hermitian 
dynamics ; 73 ! 74 possible mappings between these and our 
models deserve future studies. 
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FIG. 7: (Color online) Crossing of two groups of diabatic 
states: states |o) whose energy is time-independent and states 
1 6} whose energies are reduced with constant velocities Vb- 
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APPENDIX A: LANDAU-ZENER TRANSITION 
IN A MULTI-LEVEL SYSTEM 

We consider the multi-level Hamiltonian 

H(t) =J2s a \a)(a\+J2( £ b- v bt)\b)(b\ 

" " (Al 

+ J2(Xa b \a)(b\+X: b \b)(a\) 

a , 6 



For the experimentally important hybrid situation in 
which the bath causes both relaxation and dephasing, it 
was found that the influence of dephasing is negligible, 
unless the qubit-bath coupling is exceptionally strong. 
This robustness of the LZ transition probability under 
dephasing is quite surprising and important in applica- 
tions. Furthermore, the application of our results to ex- 
periments in superconducting circuits seems very promis- 
ing, for example for circuit QED and for the fast-pulse 
readout of phase qubits. 

In the future, it will be interesting to clarify whether 
a degeneracy of the ground state modifies the results. 
The effects of finite temperatures on the LZ tunneling 
would also deserve thorough investigations in the light 
of the exact zero-temperature results, with the recent 
work by Pokrovsky and Sun^ 3 - as an important first step. 
Other experimentally relevant issues are how nonlinear 
sweeping and finite sweeping times affect the dissipative 
transition. 



with the time-dependent diabatic energies sketched in 
Fig. [Jj We assume that all diabatic states |o), \b) are 
mutually orthogonal and that all Vb > 0. In the limit 
t — ► ±oo, the states \a), \b) become eigenstates of the 
Hamiltonian (|A1[) . The off-diagonal part of the Hamilto- 
nian is such that it only couples states of different groups 
while states within one group are uncoupled. Without 
loss of generality, we assumed in (|AL|) that the diabatic 
energies of the a-states are time-independent. If they all 
had an identical velocity v a smaller than all vj,, we could 
obtain the Hamiltonian (|A1|) by a gauge transformation 
with the time-dependent phase factor exp(— iv a t 2 /2h). 
Then vj, becomes the velocity of the 6-states with respect 
to v a . 

We choose as an initial condition that the system starts 
at t = — oo in one particular state \a). The central quan- 
tity of interest is the probability P a ~> a ' for a nonadiabatic 
Landau-Zener transition to a state \a') at t = oo. 

It is convenient to work in an interaction picture with 
respect to the diagonal part of the Hamiltonian (|Aip and, 
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thus, to apply the unitary transformation 
U (t) = J2 e ~ l£at/h \ a )( a \ +^e- ie » t/h+m » t2/2h \b)(b\. 

a b 

Then we have to deal with the interaction picture Hamil- 
tonian 

H(t) = ^e l(e "- eb)t/h+tVbt2 /2h X ab \a){b\. + U.c. (A3) 

a , 6 

In the interaction representation, the nonadiabatic tran- 
sition probability reads 



P a ^ a , = \(a'\S\a)\ 2 



(A4) 



where the S'-matrix is given by the usual time-ordered 
exponential 



k=0 



2k 



Si- 



with the (2fc)th order contribution 

Sk = 



(A5) 



dh / dt 2 --- dt 2k H{t 2k )...H{t 1 ). 

J Jti ***2fc — l 

(A6) 

Note that the interaction-picture Hamiltonian H(t) al- 
ways flips between states \a) and \b), so that the pertur- 
bation series for the transition amplitude in (|A4[) only 
consists of even powers of the interaction X ab . 

In order to compute the matrix element (a'\Sk\a), we 
will generalize and simplify our reasoning of Refs.[l8| and 
[H for this more general model. To begin with, we insert 
2k — 1 times the unit operator ^ Q \a)(a\ + J2 b an d 
obtain 



^ ^ Xab k X* k _ ibk X ak _ lbk l . . . X aibl X* bi 

i...a fc _i 6i...b fc 

f°° f°° ( i 

< dt±... dt 2k exp<-(e bl ~ e a )ti 

J-oo Jt 2k -l 1/1 

. k ft— 1 

1=2 1=1 
k 



1=1 



(A7) 



In order to cope with the time-ordering, we substitute 
the time variables ti, . . . , t 2k by the time differences Tt — 
t 2 £ — t 2 i~i and U£ = t 2 £ + i — t 2 g as sketched in Fig. [8^, i.e. 
we set 



U = t 



'=1 



'=i 



he+i = t + ^(t# + u t ) 



(A8a) 
(A8b) 

(A8c) 



(a) 



t n ui 
H 1 h 



Ujt-l 



n 



+ 



t 2 £- 2 t 2 £-\ 



i-2 r 



n\ t Ml U 2 

( b ) 1 1 h 



Ut-l 

1 1 ► 

Si S 2 S 3 Sf_i Si 

FIG. 8: Relation between the various time variables. 



which is equivalent to T( = t 2 g — t 2 g_i and ug = t 2 g + i—t 2 g. 
Note that the Jacobian of this substitution is 1. Then the 
multiple time integral in expression (|AT[) becomes 



dt I dui...du k -i / dr\...dr k 
oo Jo Jo 

k 



. k . k-i 



(A9) 



k 



=i 
i 



i=i ~ ^'=i 
Performing the i-integration, we obtain the Dirac delta 

k 



2ttHS (e ; - £a + «/3 f T£) 



(A10) 



Because all v b > and the integration interval for each 77 
is [0, 00), the sum in the argument of the delta function 
is non-negative. This has two important implications: 

First, the energies of the final and the initial state must 
fulfill the relation e a / < e a , which means that all states 
\a') with an adiabatic energy higher than that of the ini- 
tial state I a) are finally unoccupied, or 



P„ 



for e a i > e a . 



(All) 



We call this no-go theore m 56 ! 57 the "no-go-up" theorem. 
As a corollary, we find that if the system starts at t = — 00 
in the adiabatic ground state, the final state will be cither 
the initial state or one of the states \b). 

If all avoided crossings of the adiabatic energies are suf- 
ficiently narrow, the no-go-up theorem can be understood 
by semiclassical considerations: At each avoided crossing, 
the population splits up into a coherent superposition of 
two branches. If the system starts out in an a-state which 
has velocity v a = 0, it can only choose between staying 
in an a-state with constant energy or an 6-state with de- 
creasing energy. Thus transitions to diabatic states with 
higher energies are impossible. Note that the validity of 
these semiclassical arguments is limited while the no-go- 
up theorem (|A11[) is an exact statement that holds for 
any width of the avoided crossings. 
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The second implication of Eq. (|A10|) concerns the spe- 
cial transitions |o) — * \a), to which the remainder of this 
Appendix is devoted. Since in this case, the sum in the 
delta function must vanish, any contribution to the corre- 
sponding transition probability must come from the sub- 
space n = T2 = . . . = Tfc = 0. This means that the r- 
integrals in expression (|A9[) vanish unless the it-integrals 
are singular for Tg = 0. This finding can be exploited for 
the simplification of the u-integrals. In the present form, 
however, the twofold summation in the last term of the 
exponent in expression (|A9|) will complicate this task. 
Therefore it is convenient to substitute the integration 
variables u\, . . . , Ujt_i by the times 



(A12) 



which means ug — sg+i — sg as sketched in Fig. [SJd. Then 
we obtain the partially time-ordered integral 



OO /-OO /-OC poo 

dsi / ds 2 ■ ■ • / ds k / dn . . . dr k 

-00 J si J J 

. k-1 

exp \ji.£a> - £a k _ 1 )s k + ~ ^2{e ae - Sag-Jst 



1=2 



. k 



(A13) 



VbJgSg 



K £ — 1 



Setting a k = a' and ao = a, one sees that all s-integrals 
are of the form 



dse exp{i(e af - e at _ 1 + Vb t Te)se}, (A14) 



where the lower integration limit can be finite or —00. 
Evaluating this integral, one finds two types of terms: 
The first one is a principal value which is always regular 
and, thus, will not contribute to S k - A second term is 



proportional to the delta function 8(s a 



v be Te). 



This will contribute if and only if its singularity is located 
at T£ = 0, as discussed above. Therefore, we find that 
energies of all participating a-states must be identical, 



£ a - 



(A15) 



In the absence of degeneracies in the spectrum of the a- 
states finally, the important condition follows that all 
non-vanishing contributions to the perturbation series 
must fulfill the relation 



a 1 



a k 



(A16) 



This selection rule states that the only allowed processes 
are those in which the system jumps repeatedly from the 
initial state to one of the 6-states and back. Note that 



this selection rule holds only for the contributions to the 
final transition probability at time t = 00. At interme- 
diate times, other a states can be populated as well, as 
has been exemplified in a numerical study of Landau- 
Zener transitions of a qubit coupled to a single harmonic 
oscillator^ 

By use of the selection rule (|A16[) and the form (|A13|) 
of the multiple integral, we obtain for the 2fcth order term 
(a\S k \a) the expression 

E E \ X ab k \ 2 \ X ab k - t \ 2 ■ ■ -\X abl \ 2 
oi...Ofe_i bi...b k 



I dsi / ds2 ■ ■ ■ ds k / dri . . . dr k 

-00 J si J Sfc-i J 

J I2l=l Vb e TlS e /tl 



(A17) 



x e 



. k 



x exp { " \ E " £a ~ y ~ H T#)n}. 

1=1 i'=i 

A most important observation is now that the matrix 
elements X a b t no longer depend on the index ag. Any 
permutation of the integration variables sg can therefore 
be undone by proper relabelling. Thus, we can replace 
the s-integrals in the second line of Eq. (|A17|) by the 
symmetrized version 



1 f°° ■ k 

— J dsi... ds k ei 



(2irh) k 6(T 1 )...S(T k ) 

k\ v bl ■ ■ ■ v bk 

(A18) 

The remaining r-integrations can be evaluated as well, 
each of which yielding 1/2, so that finally 



Vb 



(A19) 



Inserting this into the series (|A5[) . we obtain the exact 
nonadiabatic Landau-Zener transition probability 



Pa- 



exp 



2tt x ^ 

b 



(a\X\ 



i'b 



(A20) 



where X = ^ a b \a)X a b(b\ +h.c. denotes the off-diagonal 
part of the Hamiltonian (lAljl . This general formula gen- 
eralizes our previous resultsi 18 i 38 Similar formal results, 
not applied to quantum dissipation, have been presented 
in the very recent work by Volkov and Ostrovsky.— 

Of much practical importance is the case in which 
all 6-states have the same velocity, Vb = v, so that we 
face a situation of two crossing energy bands. Owing to 
(a\X\a') = for all a, a', one finds J2 b ( a \ x \ b ) ( b \ x \ a ) = 
(a\X 2 \a). For a nondegenerate initial state \a) we end up 
withi&S 8 . 



Pa 



exp 



27r(a|A" 2 |a) 

hv 



(A21) 
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